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Abstract 

We consider the multi-channel inverse scattering problem in one-dimension in the presence of 
thresholds and bound states for a potential of finite support. Utilizing the Levin representation, 
we derive the general Marchenko integral equation for N-coupled channels and show that, unlike to 
the case of the radial inverse scattering problem, the information on the bound state energies and 
asymptotic normalization constants can be inferred from the reflection coefficient matrix alone. 
Thus, given this matrix, the Marchenko inverse scattering procedure can provide us with a unique 
multi-channel potential. The relationship to super symmetric partner potentials as well as possible 
applications are discussed. The integral equation has been implemented numerically and applied 
to several schematic examples showing the characteristic features of multi-channel systems. A 
possible application of the formalism to technological problems is briefly discussed. 
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I. INTRODUCTION 



The first attempt to extend the sinerle-channel inverse scattering problem (ISP) on the 
hne 111 12; 121, ^j. 5|| to a wider class of N x N coupled-channel potentials was made by Wadati 
and Kamijo j6| about 25 years ago. They derived a Marchenko equation associated with 
the N X N Schrodinger equation on the entire line (—00 < a: < 00). The problem has 



also been discussed in detail by Calogero and Degasperis p|. In these investigations the 
presence of threshold energies was not included and, most importantly, practical aspects of 
the implementation of the solution of the ISP were not considered. Twentvyears later, the 
coupled-channel problem was taken further with the inclusion of thresholds j8| and numerical 
solution of the corresponding Marchenko equation. 

The lack of progress can be attributed to several reasons. Firstly, the solution of the ISP 
presupposes the knowledge of the full reflection coefficient (moduli and phases) at all incident 
momenta q. This requirement can not be easily fulfilled, especially due to the missing phase 
information in standard experiments, a diffiuclty which is similar to the well known and 
longstanding phase problem in diffraction analysis ^ 10, 11 1. Secondly, a profile with an 
unlimited support leads to tough numerical questions because the highly oscillatory behavior 
of the reflection coefficient makes the numerical procedure cumbersome and unstable. A 
third difficulty in the presence of bound states is the determination of the bound state 
normalization constants for the various channels which are crucial for a unique reconstruction 
of the underlying coupled-channel potential. 

Coupled-channel inverse scattering techniques in one dimension have been considered and 
numerically implemented for specular refiection of polarized neutrons from plane surfaces of 
magnetized samples These problems contain neither thresholds nor bound states and 
can be treated by the available Marchenko equations. The formulation of coupled-channel 
inverse scattering techniques taking simultaneously into account thresholds and bound states 
would be a valuable progress, specifically for the design of semiconductor quantum devices 
and for the synthesis of quantum heterostructures with specific spectral properties [isj ]. 
In this context the spectral design of systems with specific bound states embedded in the 
continuum would be of great interest for application. 

In the present work we focus on the formulation of equations and numerical procedures 
for the solution of the coupled-channel inverse scattering problem on the line including 



thresholds. We also take into account bound states lying energetically below the thresholds 
of all channels, but exclude those embedded in the continuum of any channel. Further- 
more, we point out that the neglect of bound states in the input of the inversion procedure 
results in a supersymmetric partner potential with the same reflection properties, but com- 
pletely different spatial shape. To test the formulated inverse scattering procedure and their 
numerical implementation, model examples are considered. In these models the reflection 
coefficient has been evaluated using known profiles and employed as input in the inverse 
scattering equations to recover the original scattering profile. The agreement between the 
reconstructed and the original profile is an important measure for the quality of the inverse 
scattering procedure and a criterion for the reliability of the numerical methods. 

In Sect. II we present the formalism and give the essentials for the derivation of the gener- 
alized Marchenko equation. Technical details concerning the evaluation of the corresponding 
integrals and the treatment of bound states are outlined in the Appendix El Furthermore, 
we sketch briefiy the formalism of supersymmetric quantum mechanics for coupled- channel 
systems with thresholds and the construction of supersymmetric partner potentials. The 
feasibility of the developed formalism is demonstrated in Sect. Ill, where we reconstruct 
successfully coupled-channel potentials from simulated reflection and bound state data for 
various two-channel systems. Examples with and without thresholds and bound states are 
considered. In Sect. IV we summarize our conclusions and discuss potential applications of 
this inverse scattering technique. 

II. FORMALISM 

A. The direct problem 

The multi-channel Schrodinger equation in one dimension has the form (in appropriate 
units) 



where V(x) is assumed to be a real symmetric N x N potential matrix of finite support 
with matrix elements Vij{x), i,j = I,-- - ,N. Without loss of generality we assume for 
convenience that V(x) vanishes on the entire negative x-axis. The solution ^ is a x 
matrix whose columns are formed by the A^ linearly independent solution vectors of (P) at 




(1) 



an incident momentum k and £ is a diagonal matrix containing the threshold energies e^, 
i = 1, ■ ■ ■ ,N. We assume that the threshold energies, are arranged in increasing order 
with the lowest (in most cases the entrance channel) set to zero, ei = 0. The special case of 
no thresholds, £= 0, has already been considered by Wadati and Kamijo jgl in 1974. Here, 
we focus on the general case and allow the occurrence of non-vanishing thresholds, ej > 
forj = 2,3,...,iV. 

The free solutions are given by exp(±i/Cx), where now /C is given by 

= kH-£, 

1 being the A^- dimensional unit matrix. The matrix K. is diagonal and is defined on the 
physical sheet of the Riemann surface for the momentum variable k. This sheet has an 
(A^ — l)-fold branch cut on the real axis on the upper rim of which ki = k, and the other 
diagonal elements of K are defined as follows 



kj 



^^k'^ — tj for k > ^Jtj 

+i^/e~^^ for \k\ < ^e" (2) 
^-\^k^^j for k<-^ 



forj = 2,--- ,iV. 

Similarly to the one-channel case we introduce sets of linearly independent matrix Jost 
solutions ^±{k,x), !F±{k,x) defined by the boundary conditions 

lim exp(=i=i/Cx)^-|-(A;, x) = 1 (3) 
lim exp(±i/Cx)^-|-(A;, x) = 1. (4) 

The boundary condition !F-{k^x), for example, implies that each solution vector has at 
X — >■ — oo only in one channel an incoming wave from the left (e**^^^) and none in all other 
channels. It is therefore obvious that jF_.(/c,x) will provide the incoming component of the 
physical solution ^^(A;, x) for incidence from the left. Because the Jost solutions form a 
complete set we can write the physical solution of for left and right incidence in terms 



of J-±{k,x) and ^±{k,x), 

= ^-{k,x) + :F4k,x)nL{k), (5) 
^R{k,x) = J'4k,x)rn{k) 

= ^4k,x) + j^+ik,x)nRik). (6) 

Here, 'R'L,R{k) and T ^^^{k) are the reflection and transmission coefficient matrices for left 
(L) and right (R) incidence of the beam, respectively. 

In the absence of thresholds the members of each pair of solutions are related by replacing 
k by —k, i.e. 7-±{k,x) =^±{—k,x) (cf. ^j). In the presence of thresholds, however, this 
relation does not hold anymore for k'^ < because of the A^-fold connectivity of the /c-plane 
Considering the symmetry of the transformation (j2)) and of the Schrodinger equation 
dH) it is obvious that the identities 

^±{—k^x)=^*^{k,x) and ^±{—k,x) = ^^{k,x) (7) 

are also satisfied in the presence of thresholds. A direct consequence of Eq. ((7j) is the 
symmetry property of the reflection and transmission matrices for real k, 

nR,Li-k) = n%M and rnA-k) = r%^ik). (8) 

Introducing the matrix generalized Wronskian relation, 

WW(k,x),^k,x)] = <il^(k,x) (^^(k,x 

\ ax 



- (^-^^f^{k,x)^ ^k,x), (9) 



where \Ef and $ are solutions of the coupled- channel Schrodinger equation and T denotes 
transposition, one obtains 

4-W[^{k,x),^k,x)]=0. (10) 
ax 

This means that the Wronskian of two solutions of Eq. is constant on the entire x-axis. 
Specifically, it vanishes for two linearly dependent solutions. Evaluating the Wronskian 
ly [\E'l, \E'l] = in the limit x — oo and using the relation of \E'l(A;,x) given in Eq. Q, 
yields 

nlK = ktZl . (11) 



This relation clearly indicates the effect of the thresholds on the symmetry of the reflection 
matrix. 

B. Inverse Scattering Equation 

The determination of the potential matrix, V(x), from the knowledge of the reflection 
matrix, IZiik) (7?./j(A;)), is known as inverse scattering problem and its solution can be 
obtained via integral equations, often referred as Marchenko equations, for several quantum 
systems. Recently, an integral equation for one-dimensional coupled-channel systems with 
thresholds has been presented The derivation did not include the presence of bound 
states which are essential features of realistic quantum systems. Here, we focus on the ISP 
of coupled-channel systems in one dimension including thresholds and bound states. In the 
following the essential aspects of the derivation are outlined. 

The integral equation for the solution of the ISP associated with Eq. (HI) is most easily 
obtained via the Levin representation [21,13] of ^±{k,x), 



Inserting these expressions into Eq. leads after some algebraic manipulations to differ- 
ential equations for the transformation kernels B±{x, y), 



where we have introduced the commutator [A.,13]=yAJ3-13A.. The transformation kernels are 
related to the potential via 




(12) 



(13) 




(14) 




(15) 



and satisfy the boundary conditions 



a;,i/— >itoo 



lim B±{x, y) = . 



(16) 



Thus for X = y 




(17) 



The above partial differential equations together with the boundary conditions (fTH|l and (fTTjl 
constitute a Goursat problem 5| of a generalized nature. 

Because of the restriction placed on the potentials, the existence of solutions for the 
partial differential equation ()14|) under the above boundary conditions, can be easily shown 
as follows: A change of variables to u = |(x — y) and v = |(x + transforms p4|) into 

|2^ = V(« + .)S_ + |f.S_|; (18) 

A formal iterative solution of this partial differential equation involves repeated integrations 
over the potential. Since we imposed the condition of finite support on V this series will 
converge, due to the triangular nature of domains over which the integrations take place. 

The corresponding Marchenko integral equation is obtained by multiplying Eq. (0) from 
the right with 

1 1 

— exp(-iX:?/) — = — exp(-i/C?/)/C-^fc, (19) 
/TT dk In 

substituting the Levin representations for !F-{k,x) and ^-{k,x) on the right hand side 

and integrating over k. Restricting ourselves to systems which have only bound states at 

negative energies (fc^ = < 0, a = 1, . . . , A*";,) and assuming that no bound states are 

embedded in the continuum, we can evaluate the relevant integrals (see Appendix^. The 

result can be cast into the form 

B4x,y) + p_{x,y) + j dzB_{x, z)p_{z,y) =0. (20) 

where x > y. The input kernel of the above integral equation is given by 

P-ix,y) = p_{x,y) + p'^^\x,y) (21) 

with 

p4x,y) = dke-'^''nL{k)e-'^y , (22) 

oo 
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p^^\x,y) = -iY.^-'^-^Mae-'^-'Q-'qa, (23) 



a=l 



and Qa =K{qa). 



C. Determination of the Input Kernel 



An important step for the application of the inverse scattering procedure is the de- 
termination of the input kernel p_ (x, y) from the scattering data. Using the relation 
^hik) ='K,\{—k) one can immediately evaluate the contribution of the continuum spec- 
trum, p^{x,y). The bound state contribution, p^^\x,y) requires the knowledge of the 
bound state energies and the corresponding asymptotic normalization matrices AIq,. If 
these quantities are obtained from simulated data, assuming a given potential V(x), then 
the inversion procedure should yield the original potential. In spectral design, however, 
these quantities are not available except when one is interested to have a bound state at 
a specific energy. In such a case arbitrary values of Ada, niay result in a rather extended 
profile which is of limited interest because it can not easily be materialized. 

To ensure that the inverse scattering procedure provides us with a potential V(x) which 
vanishes for x < we must look for the required characteristics of p_ (x, y) . From the 
Levin representation of Eq. p3|) it follows that 13^{x, y) vanishes identically for x < 0. 
Because of Eq. fEUj) . the input kernel p_{x,y) must also vanish for y < x < and Eq. (j^U 
reduces for x < to the identity 



where y < x < 0. This relation is formally equivalent to the one-channel case [ij] and 
expresses the fact that the potential V(x) can be forced to vanish for x < if there is a 
series of bound states with binding energies and normalizations AIq which compensate 

Because of the peculiar spatial form of the contribution, it is obvious that those bound 
states which are energetically closest to threshold determine the behavior of the potential 
at large negative x-values. Thus for practical applications, we firstly assume to have a 
finite number, A^";,, of bound states, and restrict ourselves to those energetically closest to 
threshold. With a fixed number Nf, of bound states, it is then possible to deduce the bound 
state parameters solely from the knowledge of IZl via the nonlinear Eq. (jSHI)- We may solve 
this equation in two steps. Since the diagonal elements of B_(x, y) only depend on x + a 
one dimensional fit to the sum of exponentials will produce the g^-values and the diagonal 



P-{x,y) 



pL'{x,y) 



(24) 



oo 




a=l 



elements of Aia- The knowledge of the values is then sufficient to obtain the remaining 
non diagonal elements of the Ada- Therefore, the kernel p_{x,y) can be recovered and the 
determination of the potential can be achieved via the Marchenko equation. It should be 
noted that in the presence of thresholds, the input kernel p_(x, y) depends on the variables x 
and y separately, while in the case of no thresholds it depends on the sum x+y. Furthermore, 
the input kernel has the symmetry property p^{x,y) =p^{y,x) because of the symmetry, 
Eq. (HH), which implies 'RlIC'^ =K-^7ll. 



D. Supersymmetric Partners 

The omission of the contribution of the bound state spectrum in p_ (x, y) leads via the 
inverse scattering procedure to a potential which generates the same reflection coefficient but 
does not sustain any bound state. This feature is characterizing phase-equivalent partner 
potentials which can be obtained easily via techniques of supersymmetric (SUSY) quantum 
mechanics QES Q. The correspo„d,„g so-caUed SUSY tranrfonnat.ous are based on the 
factorization method which has been formulated to coupled-channel systems including 
thresholds ■ A compact outhne of the extension to coupled-channel systems is 

sketched in the following. 

In the factorization method the Hamiltonian of the coupled Schrodinger equation is 
written in the form 

+Voix) = At A^-£ + q\ (25) 



dx'^ 

where the factorization energy, g^, is smaller or equal to the energy qf of the lowest bound 
state. The index of V(x) indicates that it is the original potential before any transfor- 
mation. The factorization operators Aq are given in terms of the superpotential Wq{x), 

i^ = ±^ + >Vo(x), (26) 
which satisfies the nonlinear differential equation 

^Wo + Wl = Vo + £-f. (27) 

It is straightforward to show that for any solution "^olkjx) of the original Schrodinger 
equation (P), the transformation \E'o(/i;, x) leads to a solution at the same energy k'^ of the 



coupled Schrodinger equation with the potential 

Vi = Vo-2^Wo. (28) 

If we choose the factorization energy (f = q\-, one is able to eliminate the ground state from 
the spectrum of the original Hamiltonian (cf. jlfll . 12(1 12^). The reflection matrix IZi^k) 
associated with the transformed potential has changed, however, 

-Riik) = (Wo + iK.)no{k){Wo - iK,)-^ . (29) 

Here, Wq are the boundary values of Wo at x ±, respectively. 

To restore the same reflection matrix 'K.{k) a second SUSY transformation at the same 



energy q must be performed 



IS 



201 1 using the boundary values Wg = — Wi = —iQi S. 



Here «S is a diagonal matrix containing —1 in the flrst M rows and +1 in the remaining 
— M, where M < is the degeneracy of the ground state of the original system. 



III. EXAMPLES 



As a demonstration of the coupled-channel inverse scattering equations, derived above, 
we consider several schematic examples which exhibit the speciflc features of coupled-channel 
systems. For the numerical implementation of the integral equation pO|l we follow similar 
techniques as outlined in the appendix of Ref . Q . 

First we consider the case of a two-channel system without a threshold and a bound state. 
We choose the potential V(a;) to have different x-dependence in the various matrix elements. 
We take a Gaussian shape for Vii{x) with the parameters Vq = 0.1, 6 = 4, c = 1.8 and a two 
layer repulsive proflle (A^ = 2) for V22(x) with the parameters a = 0.01, Xq = 0.5 Vi = 0.08, 
Xi = 2.7, V2 = 0.05 X2 = 4.0. The off-diagonal elements Vi2{x) = V2i{x) are given by an 
Us = 3 sea-saw potential with the parameters Vq = 0.075, X£ = 1.2, Xs = 0.75, si = +1, 
and S2 = S3 = — 1. A deflnition of the shapes and the parameters is given in Appendix iBl 
With this potential we solved Eq. ((H) and evaluated Ti-L^k) up to kmax = 12. Using this 
7?.L(fc)- values we reconstructed the potential V{x) via the procedure given in Eqs. (plj) to 
()25|). The reconstructed potential matrix elements are displayed in Fig. 1 together with the 
original ones. It is seen that the reproduction of the original potentials is excellent. 

Next we consider the case where thresholds are present. Again we use an input potential 
V(x) which is chosen differently to the previous case in order to demonstrate the ability of 



the algorithm to deal with rather different situations. The potential chosen is a Gaussian 
with parameters Vq = 0.15, 6 = 9, c = 1.8 for Vii(x); a one layer repulsive profile (A^ = 1) 
with a = 0.05, xq = 1.0, Vi = 0.20, xi = 2.8 for V22(x); and Gaussian potentials with 
Vq = 0.12, 6 = 9, c = 2.2 for Vi2{x) = V2i(x). A threshold energy of e = 0.025 is assumed 
in the second channel. The reconstruction from evaluated TZ^lk) data up to k^ax = 12 are 
shown in Fig. 2 and are once again in almost perfect agreement with the original. It is 
interesting to note here that although all elements of V(x) are repulsive, the coupling of 
channels may lead to bound states. In the present example an increase of the strength of the 
coupling potential e.g. from Vq = 0.12 to Vq = 0.5 results in a bound state at Eh = —0.00053. 
This is in agreement with the observations made in Ref. j2l| . 




An example of a two-channel system with a bound state, but without threshold is shown 
in Fig. 3. The potential is chosen to be of Gaussian form with Vq = 0.15, b = 1.5, c = 2.2 
for Vii; a one-layer profile with a = 0.1, xq = 1.0, Vi = —0.1 and X2 = 3.3, for ¥22- 
The off-diagonal profiles Vi2{x) = V2i{x) are given by a sea-saw potential with Vq = 0.1, 
Xe = 1.5, Xs = 0.70, and Ug = 2, with Si = S2 = +1. The system sustains a bound state at 
Eb = -0.01558 with Mn = -0.0065, M12 = 0.0216, M21 = 0.0216, and M22 = -0.0711. As 
expected, the off-diagonal asymptotic normalization constants are, in this case, the same. 
As one can see from Fig. 3, the reproduction of the potential is, for all practical purposes, 
perfect. 

It is interesting to consider next the simultaneous existence of thresholds and bound 
states. For this we consider the previous example with a threshold e = 0.01; the rest of the 
input data remaining the same. The presence of the threshold generates a bound state at 
Eb = -0.0068 with Mn = -0.0104, Mu = 0.0400, M21 = 0.0257, and M22 = -0.1031. The 
reproduction of the potential is, once more, excellent. 

As a final example, we consider the SUSY transformations. Setting the asymptotic nor- 
malization constants equal to zero results, as expected, in supersymmetric partner profiles 
displayed in Fig. 4 which do not sustain a bound state and look considerably different from 
the original ones. 

It should be noted here that an insertion of a bound state at an arbitrary energy re- 
sults in different values for the asymptotic normalization constants. When these constants 
are evaluated via ()23|) . then the inversion procedure provides us always with the original 
interaction. 




IV. CONCLUSIONS 



We have studied the inverse scattering problem on the hne including thresholds and 
derived the corresponding integral equation of Marchenko type. Specific care was taken for 
the integrations on the physical sheet of the Riemann surface of the momentum variable 
k because of the (A^ — l)-fold branch cut as discussed in Sect. IIA. We found that these 
branch cuts do not generate additional problems as long as we deal with systems sustaining 
only true bound states. The integration contour is well defined and the integrand satisfies all 
conditions required for the application of Cauchy's integral formula. In this way the resulting 
Marchenko equation can be applied when threshold and bound states are simultaneously 
present in the system. 

Much emphasis was given to the application of the method which implies the numeri- 
cal implementation as well as the generation of the input kernel p-{x,y) from scattering 
data. Following the techniques outlined in Ref. we achieved almost perfect numerical 
reconstructions as it can be seen from the examples given in Sect. III. 

Furthermore, we considered the bound state contribution to p_ (x, y) which depends on 
the bound state energies and the so-called asymptotic normalization constants, where the 
latter are usually not accessible to experiment. We pointed out that the neglect of the bound 
state contribution results in a supersymmetric partner potential and leads to completely 
different profiles which, nevertheless, generate the same refiection coefficient matrix. Thus, 
the absence of bound state data may lead to wrong conclusions concerning the profiles. 

We proposed a method for the determination of bound state data with regard to the 
design of quantum devices with pre-designed refiection properties. Thus, a profile of finite 
range can be extracted which might be utilized by nanotechnology. Although we can not 
force the profile to be limited to a; > 0, the procedure allows to create non- vanishing profiles 
with a rather sharp edge at x = 0. 

A severe problem for applications is the necessity to provide scattering information in the 
energetically closed channels. In the examples presented here, we overcome this difficulty 
by using simulated data. However, for future applications a practically feasible and problem 
orientated method for analytical continuation has still to be developed. 

As already mentioned, the design of quantum devices via nanotechnology seems to be a 
potential field of application of these inverse scattering methods on the line (see also Ref. 



In the presence of couplings between different electronic bands, the construction of 
profiles with specific reflection and/or transmission properties via fitting methods guided by 
intuition becomes rather difficult. Here, the developed coupled-channel inverse scattering 
methods can give a reliable guide line towards the required profile for the device. In this 
context the inclusion of bound states in the continuum in the inversion scattering procedure 
would be of great interest because slight changes in the retrieved profiles may result in sharp 
resonances which could be of interest for applications. Work is in progress to handle this 
important question from the inverse scattering point of view. 

The presented formalism can be applied to any technological problem where one wants 
to retrieve or design the profile from the spectral information. Nanostructure devices and 
waveguides are obvious examples where these procedures can be of great interest. 
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APPENDIX A: DERIVATION OF THE MARCHENKO EQUATION 

The formulation of the solution of the ISP by means of an integral equation, frequently 
referred to as the Marchenko equation, is most easily achieved via the Levin representation 
of the Jost solution ^_{k,x) (cf. Eq. (UHl))- In a first step we multiply Eq. (0) from the 
right by 

1 1 

— exp(-iX:y)— = — exp(-i/C?/)X:"^fc, (Al) 

substitute the Levin representations for ^^{k,x) and !F-{k^x) on the right hand side and 
integrate over k. Writing T l = l+r(A;) and reordering the terms leads to the relationship 

Xi + X2 = p_ix,y) (A2) 

/X rx 
dzB-{x, z)p_{z,y) / dzB-{x, z)T-l{z,y) 
-oo J —oo 



with 



/ dkkJ^+{k,x)T{k)e-''^yK- 



J — oo 



— r dkk[J^+{k,x)-e-'^'']e 
27r J^^ 





dkke-''^''nLik)e-''^yK-^ 



dkke-'^^^'-y^K-^ . 



In the presence of thresholds these integrals must be handled consistently and with care 
because of the (A^ — l)-fold branch cut of the physical sheet of the Riemann surface for 
the momentum variable k. Because of the identity ()A1|) the j-th columns of the matrices 
Ii, X2, p, contain effectively integrals over the channel wave numbers kj. It is therefore 
important to consider the consequences of the mapping of k to kj, as given in Eq. ((21). As 
displayed in Fig. 5 the integral over the real axis of k is transformed into a contour extending 
also to the positive imaginary axis of kj,j = 2, . . . , N. 

To show the effect of this deformation of the contour we consider the matrix T-C{x,y) 
which is diagonal with elements 



with j = 1, ■ ■ ■ , N , kii = k. Performing a transformation from k to kj in the jj-element 
leads to a splitting of the integral 



Because of the symmetry of the integrand the second term of Eq. flA4|) vanishes, while the 
first term yields the (5-function. Hence we obtain 1-C{x,y) = 6{x — y)l. 

Making use of the Levin representation of ^-.{x,y) we can apply the expression for 
7i{x,y) to evaluate X2{x,y). Restricting ourselves to x > y, the matrix X2{x,y) vanishes 




(A3) 




(A4) 



and Eq. ()A2|) takes the simple form 



Ii{x,y) = p^{x,y) + B^{x,y)e{x -y) 




(A5) 



where Q{x) is the Heavyside function. 

In Eq. ()A5|) we are only left with the evaluation of the integral Ii{x,y). To evaluate 
Ii{x, y) we close the integral in the upper-half plane of k in order to apply Cauchy's theorem. 
To do so we have to consider the integrand with regard to its analyticity. It is well known 
from scattering theory ^3] that the occurrence of bound states leads to a pole of T{k) 
and 'K-{k) on the positive imaginary axis of k. This is best seen from Eq. (0) where the 
requirement of a normalizable solution at /c = g = m, /t > requires the relationship 

:F+(m,x) =:F_(m,a;) lim 7?.l(A;)T^^(A;) . (A6) 

k^iK 



Employing a method similar to the one used in 
T lik) at qa = ina are simple, i.e. 



it can be shown that the poles of 



TLik) = t:^^o. + ■ ■ • , (A7) 

where A/'a is the residue of at k = qa and therefore also of F^. 

In the following we restrict ourselves to systems which have only true bound states at 
/c^ = < 0, a = 1, . . . , A*";,. Thus the Jost solution is well defined in the upper half-plane of 
each channel wave number kj and the integrand does not exhibit irregular features along the 
contour of integration. In order to evaluate the integral over the upper half circle we must 
consider the behavior of the jF(/c,x) for \k\ oo. From Eq. Q it immediately follows that 



Im kj = Im k 



1 + -TTTTT + O ' ^ 



Hence, the dominant term for the vanishing of the Jost solution in the upper half plane of k 
for ^ oo is independent of the channel and leads to a vanishing of the integral over the 
half circle for x > y. Thus the integral Xi{x,y) is simply given by the sum over all bound 
state poles 

Xi = i ^ jF+(gQ,, x)Afa exp(— i Qay)Q{qar'qa, (A8) 

a=l 

where we assume that the potential V{x) sustains A^^^ bound states. 

Using the specific relationship, Eq. ()A6j) . between the Jost functions JF+ and at the 
bound states, k = qa, and entering the Levin representation for x) leads finally to 

X,ix,y) = -p^^\x,y)- r dzB.ix,z)p^^\z,y) (A9) 




(AlO) 



a=l 



where Aia are the residues of TZiik) at qa- Introducing the total input kernel 



P {x, y) = p (x, y) + pL' (x, y) 



(All) 



leads to the Marchenko equation ((201) ■ 

APPENDIX B: DEFINITION OF POTENTIAL PROFILES 

In the examples we use potentials, whose elements are given by an x-dependence of 
closed form. In the following we summarize the definitions of the potential forms and their 
parameters. 

1. Gaussian Profile 



1 + exp((xj — x)/a) J 

where N is the number of layers, a is the diffuseness parameter, and Vi the strength of the 
layer with width Xi — Xi_i with xq being the left boundary of the multilayer profile. 



V{x) = 1/oexp(-6(x-c)2) 



(Bl) 



2. Multilayer Profile 




1 




3. Multiple Sea-Saw Potential 



The sea-saw potential is defined by 



2 H 



X < Xq 

(-^ -^m— l) -^m— 1 ^ ^ -^m 





(B3) 



where Xm — + mxi and Xm — Xg -\- {m — \)x(_^ m — 1, . . . , n^. Here, Ug is the number 
of sea-saw peaks of length Xi^V^is the strength of the potential, and Sm is the sign of the 
corresponding term. The nonvanishing part of the profile is shifted by an amount Xg from 
the origin. Thus it extends between xq and Xn^. 
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FIGURE CAPTIONS 



Figure 1 

Reconstructed potential (dashed line) of a 2-cliannel system without thresholds or bound 
states. The details of the original potential (sohd line) are given in the text. 

Figure 2 

Reconstructed potential (sohd hne) of a coupled 2-channel system with a threshold €2 — 
0.025 in the 22-channel. The details of the original potential (dashed hne) are given in the 
text. 

Figure 3 

Reconstructed potential (dashed line) of a two-channel system with a bound state but with- 
out thresholds. The original potential (solid line) is shown for comparison. 

Figure 4 

Input potentials (sohd line) and the SUSY transformation potentials (dashed hne) . 
Figure 5 

Mapping of k to kj and the corresponding deformation of the contour for integration. 
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